ADHESIVE CONTACT DELAMINATING AT MIXED MODE, ITS THERMODYNAMICS 
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ABSTRACT. An adhesive unilateral contact between visco-elastic heat-conductive bodies in linear Kelvin- 
Voigt rheology is scrutinised. The flow-rule for debonding the adhesive is considered rate independent, uni- 
directional, and non-associative due to dependence on the mixity of modes of delamination, namely of Mode 
I (opening) and of Mode II (shearing). Such mode-mixity dependence of delamination is a very pronounced 
(and experimentally confirmed) phenomenon typically considered in engineering models. An anisothermal, 
thermodynamically consistent model is derived, considering a heat-conductive viscoelastic material and the 
coupling via thermal expansion and adhesion-depending heat transition through the contact surface. We prove 
the existence of weak solutions by passing to the limit in a carefully designed semi-implicit time-discretization 
scheme. 
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1. Introduction 

Nonlinear contact mechanics is an important part of mechanical engineering and receives still grow- 
ing attention due to its numerous applications. We focus on the modelling and analysis of an inelastic 
process called delamination (sometimes also referred to as debonding), between two elastic bodies, glued 
together along a prescribed delamination surface. "Microscopically" speaking, some macromolecules in 
the adhesive may break upon loading and we assume that they can never be glued back, i.e., no "healing" is 
possible. This makes the process unidirectional, viz. irreversible. On the glued surface, we consider the de- 
lamination process as rate-independent and, in the bulk, we also take into account rate-dependent inertial, 
viscous, and heat-transfer effects. The ultimate phenomenon counted in engineering modelling (and so far 
mostly ignored in the mathematical literature), is the dependence of this process on the modes under which 
it proceeds. Indeed, Mode I (=opening) usually dissipates much less energy than Mode II (=shearing). The 
difference may be tens or even hundreds of percents, cf. [ 1 20 .2 1 .43 1 . Moreover, the delamination process 
rarely follows such pure modes: in general, the mixed mode is favoured. Microscopically, this difference 
is explained either by some roughness of the glued surface (to be overcome in Mode II but not in Mode 
I, cf. IfTTT ) or by some plastification, either in the adhesive or in a narrow strip around the delaminating 
surface, before the delamination itself is triggered (and, as usual in plasticity involving trace-free plastic 
strain, again only shearing in Mode II can trigger this plastification, not opening in Mode I, cf. II20II471 ). 

In this article, we focus on a standard engineering model which, however, was never rigorously anal- 
ysed so far. Even for the isothermal variant of this model, existence of solutions has not yet been proved, 
although computational simulations are routinely launched and successfully used in applications. In par- 
ticular, we extend the analysis of lfT8l[3T1 to the mixity-sensitive case. For other results on models for 
rate-independent adhesive contact, we refer to 12311401 . The analysis of models featuring a rate -dependent 
evolution for the delamination variable was carried out in [4 - 8 28 29], cf. also the monographs 1 10 42) for 
further references. 

The initial-boundary-value problem we are going to analyse is written down in its classical formulation 
in Section |2] In the following lines, we just highlight the main features of the model, in particular focus- 
ing on the mixity of delamination modes. We confine ourselves to small strains and, just for the sake of 
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notational simplicity, we restrict the analysis to the case of two (instead of several) bodies + and fi_ 
glued together along the contact surface T c . The material in the bulk is taken to be heat conductive, and 
thus the system is completed by the nonlinear heat equation in a thermodynamically consistent way. The 
contact surface is considered infinitesimally thin, so that the thermal capacity of the adhesive is naturally 
neglected. The coupling of the mechanical and thermal effects thus results from thermal expansion, dissipa- 
tive/adiabatic heat production/consumption (depending, in particular, on the mode mixity on T c ), and here 
also from the possible dependence of the heat transfer through the contact surface T c on the delamination 
itself, and on the possible slot between the bodies if the contact is debonded. 

We consider an elastic response of the adhesive, and then one speaks about adhesive contact (in contrast 
to brittle contact where a mixity-dependent model seems to be particularly difficult to analyse). The elas- 
tic response in the adhesive will be assumed linear, determined by the (positive-definite) matrix of elastic 
moduli A. At a current time, the "volume fraction" of debonded molecular links will be "macroscopically" 
described by the scalar delamination parameter z : T c — >• [0, 1], which can be referred to the modelling 
approach by M. Fremond, see II13II141 . The state z(x) = 1 means that the adhesive is still 100% unde- 
stroyed and thus fully effective, while the intermediate state < z(x) < 1 means that there are some 
molecular links which have been broken but the remaining ones are effective, and eventually z(x) = 
means that the surface is already completely debonded at a; G r c . In some simplification, based on the 
Griffith concept lfl6l . it is assumed that a specific phenomenologically prescribed energy a (in J/m 2 , in 
3-dimensional situations) is needed to break the macromolecular structure of the adhesive, independently 
of the rate of this process. Thus, delamination is a rate-independent and activated phenomenon, governed 
by the maximum dissipation principle, and we shall accordingly consider a rate-independent flow rule for 
z. 

Let us now comment on the main new feature of the model presently analysed, i.e. its mixity-sensitivity. 
An immediate reflection of the standard engineering approach as e.g. in 01714411451 is to make the activation 
energy a = <x(V>g) depend on the so-called mode-mixity angle tpQ, For instance, if v c = (0, 0, 1) at some 
x G r c (with v c the unit normal to T c , oriented from il + to £!_), and A = diag(K n , ft t , « t ), the mode- 
mixity angle is defined as 



where |it] t and [«]]„ stand for the tangential and the normal traction. They give the decomposition of the 
jump of displacement across the boundary T c as [itj = [ujn^c + Jujt, with [uj n = Jit] • v c . In fact, to 
satisfy natural analytical requirements (viz. the continuity of the function cii below, cf. ( 11.31 )), one should 
rather consider a suitable regularization of fll.lt to avoid discontinuity at 0. For example, it is sufficient to 
take 



The coefficient K t is often smaller than K n , and a typical phenomenological form of a( ) used in engineering 



In ( 11. 2b . a\ = a(0) is the activation threshold for delamination mode I and A is the so-called delamination- 
mode-sensitivity parameter. Note that a moderately strong delamination-mode sensitivity occurs when the 
ratio an/ai is about 5-10 (where an = a(90°) is the activation threshold for the pure delamination mode 
II). Then, one has A about 0.2-0.3; cf. |43l 

In the thermodynamical context, the energy a needed for delamination is dissipated by the system in 
two ways: one part a\ is spent to the chaotic vibration of the atomic lattice of both sides of the delam- 
inating surface T c , which leads "macroscopically" to heat production (cf. also l40l Remark 4.2]), while 
another part ciq is spent to create a new delaminated surface (or, "microscopically" speaking, to break the 
macromolecules of the adhesive). Thus a = oo + oi. Consistently with the dissipation, Mode II also 
heats up considerably more than Mode I, as experimentally documented in ||30l . In view of ( 11.21 ). an option 
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suggested already in l38l is 

a (M):=ai, oi(H) := aitan 2 ((l-A)7/> G (H )), (1.3) 

meaning that plain delamination does not trigger heat production at all, and only the additional dissipation 
related with Mode II contributes to the heat production on the delaminating surface. We summarize the 
features of these particular modes in Table Q] 

The mathematical difficulties attached to the analysis of the PDE system for the present mixity-sensitive 
delamination model arise both from the proper thermodynamical coupling, and from hosting an inelastic 
rate-independent process on T c . Models combining thermal and rate-independent effects have already been 
successfully analysed in 13511 for inelastic processes in the bulk, and in iTJTl for surface delamination. The 
essential ingredient for the analysis is the satisfaction of the energy balance. In this direction, the concept of 
energetic solutions to rate-independent systems recently developed in [22 24-26 ] and adapted to systems 
with inertia and viscosity in 11341 appears truly essential. Here, additional difficulties are related with the 
mixity-dependence of the dissipation, which makes it non-associative, in contrast to the mixity-insensitive 
case and to another model recently devised and analysed in an isothermal case in [38, 39]. This analytical 
feature has led us to resort to a higher-order gradient in the momentum equation via the concept of the so- 
called hyperstresses, already justified and used in the theoretical-mechanical literature, cf. e.g. 115112711461 . 
Such a regularization brings various inevitable technicalities into the classical formulation of the problem, 
cf. (O and ( 1331 ) later on. 

The main result of this paper is the existence of solutions to the initial-boundary value problem as- 
sociated with the mixity-sensitive model under investigation. The proof is performed by passing to the 
limit in a suitably devised semi-implicit discretization scheme, cf. 15.41 . Let us mention that such a kind 
of scheme (already applied in |fl9l Sect. 3.1] for a special dynamic isothermal fracture problem) leads to 
considerable analytical simplifications, in comparison with the fully implicit scheme used in J3T|. In the 
existence proof we shall distinguish the dynamic case, involving inertial terms in the momentum equation, 
and the quasistatic one, where inertia is neglected. In the latter situation, we will be able to tackle a fairly 
general contact conditions for the displacement variable u, in particular including (a generalization of) the 
Signorini frictionless contact law. For further explanations and comments, we refer to Remark l43l 

The plan of the paper is as follows: in Section|2]we formulate the initial-boundary-value problem in its 
classical formulation. We also briefly sketch its derivation, referring to OTl for more details. Hence, in Sec. 
after introducing a suitable transformation for the heat equation, we advance a suitably devised weak 
formulation of our PDE system, and comment on its relation to the classical formulation of Sec. [2] We state 
our main existence result in Sec. [4] and set up the approximation via a semi-implicit time-discretization 
scheme in Sec. [5] For the discrete solutions, suitable a priori estimates are obtained, which allow us to pass 
to the limit in the time-discrete approximation, and conclude the existence of solutions in Sec. [6] 

2. The model and its derivation 

Hereafter, we suppose that the elastic body occupies a bounded Lipschitz domain Q C M 3 . We assume 
that 

n = fi + ur c ui]-, 

with f2 + and il_ disjoint Lipschitz subdomains and T c their common boundary, which represents a pre- 
scribed surface where delamination may occur. We denote by v the outward unit normal to dil, and by v c 
the unit normal to T c , which we consider oriented from f2 + to fi_. Moreover, given v £ W 1,2 (£l\T c ), v + 
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(respectively, v ) signifies the restriction of v to fl + (to resp.). We further suppose that the boundary 
of il splits as 

on = r D u r N , 

with r D and T N open subsets in the relative topology of dft, disjoint one from each other and each of them 
with a smooth (one-dimensional) boundary. Considering T > a fixed time horizon, we set 

Q := (0, T) x 0, S := (0, T) x 0Q, S c := (0, T) x T c , E D := (0, T) x T D , E N := (0, T) x T N . 

For readers' convenience, let us summarize the basic notation used in what follows: 
u : f2\r c — > R 3 displacement, 
8 : f2\r c —> (0, +oo) absolute temperature, 
z : r c — > [0, 1] delamination variable, 
e = e(u) = iVit T + |Vu small strain tensor, 
[it] = u + |r c — u~ |r c jump of it across r c , 
a stress tensor, 

$ rescaled temperature (enthalpy), 
C G R elasticity constants, 
B G R 3 viscosity constants, 
HI G R 3 elasticity constants for hyperstress, 
Gel 3 viscosity constants for hyperstress, 
K=K(e, #)gR 3x3 heat-conduction coefficients, 
E G R 3x3 thermal-expansion coefficients, 
B := CE, 

The state is formed by the triple (u, 8, z). We use Kelvin-Voigt's rheology and thermal expansion. As a 
further contribution to the stress a : (0, T) x 17 — > R 3x3 , we also consider the hyperstress f) : (0,T)x!]4 
R 3x3x3 , which accounts for "capillarity-like" effects. Mimicking Kelvin-Voigt's rheology, we incorporate 
in the hyperstress contribution to a the corresponding dissipation mechanisms. Hence we assume the stress 
a : (0, T) x Q -> R 3x3 in the form: 

a = a(u, u, 6) := Be{u) + C(e(u)-E0) - div( HVe(u) + GVe(il) ), (2.1) 

viscous elastic =: f) 

stress stress hyperstress 

The G-term will ensure the acceleration u to be in duality with the velocity ii, which is needed if q > is 
considered. Furthermore, we shall denote by T = T(u, v, 8) the traction stress on some two-dimensional 
surface T (later, we shall take either V = T c or T = T N ), i.e. 

T(u, u, 8) := a(u, u, 8) \ v , (2.2) 

where of course we take as v the unit normal v c to r c , if T = T c . 

We address a generalization of the standard frictionless Signorini conditions on T c for the displacement 
u. This is rendered through a closed, convex cone K(x) C R 3 , possibly depending on x G T c . In terms of 
the multivalued, cone-valued mapping K : T c =4 R 3 , the boundary conditions on T c can be given in the 
complementarity form as 

M h o, ] 

T{u,u,8) t 0, > onr c . (2.3) 
T(u,u,8)-lul = J 

In ( 12.31 l. >z is the ordering induced by the mapping K : T c =S R 3 , in the sense that, for Vi, v% : T c — > R 3 , 

v\ >z v-x if and only if V\{x)— V2(x) £ K(x) for a. a. x G r c . (2.4) 

Likewise, >^ is the dual ordering induced by the negative polar cone to K, viz. for Ci , C2 : T c —> M 3 , 

Ci t C2 if and only if (i(x)-v > (2(x)-v for all v G K(x), for a.a. x G T c . 



A G E 3x3 elastic coefficients of the adhesive, 

q > mass density, 

c v = c v (#) heat capacity, 

eto = a o([[^]) energy stored on Fc, 

ai = ai([u]) energy dissipated on Tc, 

7/=rj([u] , z) heat-transfer coefficient on Fc, 

F : Q ->■ R 3 applied bulk force, 

wd prescribed boundary displacement, 

/ : En — > R 3 applied traction, 

G : Q — > R bulk heat source, 

g : E — !> R external heat flux, 

?/>'s (bulk and surface) free energies, 

£'s (pseudo)potentials of dissipative forces, 

f 's (bulk and surface) rates of dissipation. 
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Possible choices for the cone-valued mapping K : T c =4 R 3 are 

K(x) = {v e M 3 ; v ■ v c (x) > 0} for a.a. x £ T c , or (2.5a) 

K{x) = {vE M 3 ; v ■ v c {x) = 0} for a.a. x £ T c . (2.5b) 



In the case (I2.5ab . (12.3b reduce to the standard unilateral frictionless Signorini contact conditions for the 
normal displacement. The case (12.5bl i prescribes the normal jump of the displacement, variable at x € T c , 
to zero. Thus, it only allows for a tangential slip along T c . This may be a relevant model under high 
pressure, when no cavity of T c can be expected anyhow. Such a situation occurs, e.g., on lithospheric 
faults deep under the earth surface. Note that, in (I2.5bb . K(x) is a linear manifold for a.a. x <G T c . Indeed, 
later on we shall have to assume the latter property in the case g > 0. 

Classical formulation of the adhesive contact problem. Beside the force equilibrium coupled with the heat 
equation inside Q\E C and supplemented with standard boundary conditions, we have two complementarity 
problems on E c . Altogether, we have the boundary-value problem 



gu — cliv^lLSe^iiJ + <L.e(u) — lav — div r)J = £ , i) = Inly e(u)+hr\/ e(u) 


in Q\Zjc, 


(z.oa) 


c v (0)0 - div(K(e(u), 0)V0) = Be(u):e(u) - 9M:e(u) + GVe(u):Ve(u) 


in Q\E C , 


(2.6b) 


u = 


on E D , 


(2.6c) 


T(u,u,6)-div a (t)-v) = f 


on E N , 


(2.6d) 


(K(e(u),0)V6)u = g 


on S, 


(2.6e) 


§:(i>®v) = 


on E D U E N , 


(2.6f) 


[Be(tt) + Ce(u) -M6- div h] v c - div s ( [h] u c ) = 


on E c , 


(2.6g) 


\) + :{v c ® v c ) = f)":(^ c ® u c ) = 


on E c , 


(2.6h) 


[n] h0 


on E c , 


(2.6i) 


T(u,u,0) -div s (f) ■ v) + zk{v\~-za' {{v\) > 


on E c , 


(2.6j) 


(T(u,u,6)-div s (t) ■ is)+zA[u]-za' ([u])) ■ [u] =0 


on E c , 


(2.6k) 


'z < 


on E c , 


(2.61) 


d<oo([u]) + oi([«]) 


on E c , 


(2.6m) 


z(d-a ([u])- ai ([u})) =0 


on E c , 


(2.6n) 


dedl m (z) + $A[u}-[u] 


on E c , 


(2.6o) 


\ (K(e(u),9)V6\+ c + K(e(«), 0)V0|F C ) ■ v c + r,( [«] , z) [0] = 


on E c , 


(2.6p) 


[K(e(u),0)V0] ■v c = -ai([u])z 


on E c , 


(2.6q) 



where : denotes the tensorial product involving summation over 3 indices, and we have used the notation 
f) + = fy|n+ an d f) = t)|n_- In d2.6dl i, we denoted by div s the two-dimensional "surface divergence", 
defined by div s := tr(V s ), where tr is the trace operator (of a 2 x 2 matrix), and V s denotes the tangential 
derivative, defined by V S V = VV- f^ZA As to the involved tensorial symbols, 

K = K(e,0) is a 2nd-order positive definite tensor, (2.7a) 
i.e. a 3 x 3-matrix, while 

C, D : M 3 *^ -> M 3 * 3 are 4th-order positive definite and symmetric tensors, C potential, (2.7b) 

in the sense that both Ce(u) and Oe(it) are R 3 * n 3 -valued, and the operator divCe(it) has a potential. Anal- 
ogously, for the higher-order terms we suppose that 

H,G : R 3x3x3 -> r3x3x3 are 6th-order posit, def. and symm. tensors, H potential, (2.7c) 

in the sense that both divEIVe(it) and divGVe(u) are R 3 ^-valued, and the operator div 2 HVe(u) has a 
potential. 
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The complementarity problem ( |2.6it -( |2~6kl i describes general, possibly unilateral (depending on the 
choice of the mapping K : T c =4 R 3 ) contact. Indeed, for \) = and z — (i.e., no hyperstress 
contribution, and complete delamination), ( I2.6il )-( l2~6kl ) reduce to relations ( 12.31 ). which in turn generalize 
the Signorini conditions. For later reference, we point out that the complementarity conditions ( 12.6ib -( l2~6kT i 
are equivalent to the subdifferential inclusion 

dI K (lul)+T(u,u,6) + zAlu} -za^luj) 30 on E c , (2.8) 

which features the indicator functional Ik ■ L 2 (T C ;M. 3 ) — > [0, +oo] associated with the multivalued 
K : T c =t M 3 , viz. 

I K (v)= [ I K(x) {v(x))dS foralH>eL 2 (r c ;R 3 ), (2.9) 

and its subdifferential (in the sense of convex analysis) 81k ■ L 2 (T C ; R 3 ) =4 L 2 (T C ; R 3 ). 

In turn, adhesive contact results from the complementarity conditions d2.6U - d2~6ol l, which can be refor- 
mulated as the flow rule 

fl/(-oo,o](i) + W[o,i](«) + 3A[u].[«] -a (H)-ai([u])9 in E c , (2.10) 

with the indicator functions /(_oo,oi) ^[0,1] : R — > [0, +oo] and their (convex analysis) subdifferentials 

0J(-oo,o], dI [0A] :R^R. 

Some comments on the derivation of the model. In OTl . a thorough derivation of the analogue of system 
(12.6b . in the case when the dependence on the delamination modes is included in the model, has been 
developed. Therefore, we refer the reader to OTl Sec. 3] for all details, and in the next lines we just 
highlight the main differences between the system considered in BP and the present ( 12.6b . Namely, here 

the free energy is enhanced by the higher-order term iHVe(u):Ve(tt), the dissipation energy in the bulk 

is enhanced by GVe(«):Ve(u), and on the delamination surface T c we have ao and ai depending on [it] . 

More specifically, we consider the free energy, the dissipation rate, and the (pseudo)potential of dissi- 
pative forces in the bulk given by 

i i rfl 

^ buIk (e, Ve, 0) = -C(e-E0):(e-E0) + -HVe:Ve - yB:E - tp (6), (2.11a) 

£ bulk (e, Ve) = 2( 2 (e, Ve) with C 2 (e, Ve) := ^ID>e:e + -GVe:Ve, (2.11b) 

where ipo : (0, +oo) — >■ R is a strictly convex function. The free energy and the dissipation rate on the 
contact surface are 

^ (z^Avv-aoiv)) if^Oand 0<z<l, ^ 

I +oo otherwise, 

e ^ {v ,z )=Cl{v , i):= hM\*\ ifz<0, (212b) 
I +oo otherwise. 

The overall free energy and the (pseudo)potential of dissipative forces are then 

*(u,j?,0)= / 7/> bulk (e(u),Ve(u),0)dx + [ i/; surf ( {uj , z) dS , (2.13) 
Ja\r c Jr c 

H(u; «,*)=/ e ulk {e(u),Ve(u))dx+ [ e url ([uj,i)dS, (2.14) 
Jn\r c Jr c 

respectively. Considering the specific kinetic energy ^g \v\ 2 (with g > the mass density), for all v E 
L 2 (f2) we define the overall kinetic energy Tki n and the external mechanical load L by 

T^{v):= \ t g\v\ 2 dx and (L(t),u) = [ F(t)-udx + [ f-udS. (2.15) 

The mechanical part of system ( 12.6b . i.e. equations d2.6b .c,d,f-o). is then just the classical formulation of 
the abstract evolutionary system 

dT kia (u(t)) + (ii4) 3(u(t); e(«(t)), z{t)) + d {u ^(u(t), z(t), 0(t)) 3 L(t) for t e (0, T), (2.16) 
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where d denotes the (convex analysis) subdifferential of the functionals S and ty, w.r.t. suitable topologies 
which we do not specify. The remaining equations in ( 12.6b yield the heat-transfer problem, i.e. (12.6b .e.p.q). 
Its derivation proceeds standardly from postulating the entropy s by the so-called Gibbs' relation s = 
-%(u, z, 9), viz. (s, 9) = z, 9; 9) for all 9, where %(u, z, 0; 9) is the directional derivative of * 

at (it, z, 9) in the direction 9. This yields the entropy in the bulk as 

s = s(9,e) = ^-(e(u),6) =M:e(u) + i>' (9). (2.17) 

Further, we use the so-called entropy equation 

9s = <e bulk (e(u)) -divj. (2.18) 

Substituting s = M:e(u) — i/j'q(9)6, cf. ( 12.17b . into the entropy equation ( 12.181 ) yields the heat equation 

- div j = 2C 2 (e(tt)) - 9M:e{u) with c v ((9) = 9^(9). (2.19) 



Hence, assuming the constitutive relation j := — K(e(it), 0)V# for the heat flux, i.e. Fourier's law in an 
anisotropic medium, one obtains the heat equation in the form (I2.6bb - 

Similar, but simpler thermodynamics can also be seen on the contact boundary by involving ■0 surf and 
£ su . Since d2.12a| ) is independent of temperature, the "boundary entropy"=— ^V' surf is simply zero, and 
the corresponding entropy equation reduces to = £ surf (z) — [j] (as an analog of ( 12.181 l). which then 
results in (2.6qi. Incorporating the analog of the phenomenological Fourier's law, we arrive at (2.6p I. 

We emphasize that the model is thermodynamically consistent, in the sense that it conserves the total 
energy, i.e. here 

||u| 2 + \ce(u):e{u) + ijHV e(u)':V e(u) + h{9) dx 



./n\r c 2 



kinetic, elastic, and thermal energies 

+ f ^A[tt]-[«] -a ([u])zdS = [ F-udx + [ gdS+ [ f-udS , (2.20) 

Jr c 2 Jn t Jr Jr N 

mechanical energy power of bulk power of surface heat 

in the adhesive mechanical load an{ j mechanical load 

and it satisfies the Clausius-Duhem's entropy inequality: 

as well as non-negativity of the temperature under suitable natural conditions, cf. Theorem l4.2l 



3. Enthalpy transformation and energetic solution 

In what follows we are going to tackle a reformulation of the PDE system ( 12.6b , in which we replace the 
heat equation ( I2.6bb with an enthalpy equation. Namely, we switch from the absolute temperature 9, to the 
enthalpy §, defined via the so-called enthalpy transformation, viz. 

i? = h{9) := [ c v {r)dr. (3.1) 
Jo 

Thus, h is a primitive function of c v , normalized in such a way that h(0) = 0. We will assume in ( 14. lab 
below that c v is strictly positive, hence h is strictly increasing. Thus, we are entitled to define 

frV) if^>o, , . = x( e ,ew) 

U[ h \0 if t9 < 0, Me,v). Cv(0(i?)) , V-V 

where /i" 1 denotes the inverse function to h. 
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With transformation A3. U and the related (13. 2\ , also taking into account the subdifferential formulations 
and (I2.101 i. the PDE system (12.6b turns into 

giZ - div(Be(«)+Ce(«)-B6(i9)-div h) = F, h = HVe(u) + GVe(ii) 

j? - div(3C(e(u), #)W) = Be(u):e(tt) - e(i?)B:e(u) + GVe(i):Ve(«) 
u = 

(OC(e(u),i?)Vi?)t/ = fl 
T(u,u,i?)-div s (h •!/) = / 

f):(f (8 f) = 

pJe(u) + Ce(«) - B6(i?) - div fj]j/ c - div B ([lj]i/c) = 
h+:(^ c ® i/ c ) = t)~:(v c <S)v c )=0 

dI K (luj)+7(u,u,{)) - div s (f) • v) + zA[u] - za'odu}) 3 
W(-oo,o](«) + dl [0 ,i](z) + iAH-H - «o(M) - ai(M) 3 
^(X(e(u),^\+ c +%(e(u),^^ c )-Pc + v(H,z)lQm =0 
[3C(e(«),T?)Vi?] -^ c = -ai(H)^ 




> on S r 



(3.3a) 
(3.3b) 
(3.3c) 
(3.3d) 

(3.3e) 



with 



(3.3f) 



T(u,v,i?) :=T(u,«,6(^)) = [Be(w) + Ce(u) - B9(tf) - div(HVe(w) + GVe(w))] \ T v 
where again we take as v the unit normal to r c if T = T N , and v = v c if T = T c - 

Remark 3.1. The reformulation of the heat equation ( 12.6bl i. viz. the second of ( 13.3al i. shows the advantage 
of the enthalpy transformation ( 13. 11 1. Indeed, by means of ( 13.11 ). the nonlinear term c v (6)8 in ( I2.6bl i has been 
replaced by the linear contribution d. This makes the second of d3.3ab amenable to the time-discretization 
procedure developed in Section [5] Such a procedure would be more troublesome, if directly implemented 
on the heat equation d2.6bl ). also taking into account the growth conditions on the heat capacity function 
c v , which we shall impose in ( 14. lbb (for analogous assumptions, see H31H33l[35l and iflOl Sect. 5. 4. 2] for 
contact problems in thermo-viscoelasticity). 



Functional setup. Throughout the paper, we shall extensively exploit that, in our 3-dimensional case, 
W 2 ' 2 (Q) C W 1,p (fi) continuously for 1 < p < 6, and 



u i-> u\ r : VK 2 ' 2 (0) -> C(T) compactly, 



(3.4) 



with r = dil, or r = r c , or r = T N . Moreover, for 7 £ [2, 00) we will adopt the notation 

I 3 ) := {v £ iy 2 ' 7 (fi\r c ; R 3 ) : v = on T D } , 



and denote by (■, ■) the duality pairing between the spaces W^(Q\T C ; R 3 )* and W 2 ^(Q\T C ; R 3 ). Fur- 
thermore, in the case K(x) is a linear subspace of R 3 for almost all x £ T c , we shall use the notation 



K 



: (fi\r c ;M 3 ) := {v£W^(n\T c 



[v(x)J £ K(x) for a.a. x £ T c }. 



(3.5) 



We will work with the space of measures M(f2) := C(fi)*. Finally, let X be a (separable) Banach space: 
we denote by M([0,T];X), B w *([0, T}\ X), and BVQO, T]; X), respectively, the Banach spaces of the 
measures on [0, T] with values in X, of the functions from [0, T] with values in X that are bounded and 
weakly* measurable (if X has a predual), and of the functions that have bounded variation on [0, T] . Notice 
that the functions in B w *([0, T];X) and BV([0, T];X) are defined everywhere on [0, T]. 
Loading qualification. Hereafter, the external mechanical and thermal loading F, /, and g will be quali- 
fied by 

F £ L\0,T;L 2 {Q;R 3 )); (3.6a) 
f £W^(0,T;L^ 3 (r N ;R 3 )), (3.6b) 



g£L 1 (Y>), .g>0a.e. in E 



(3.6c) 
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Initial conditions qualification. As for the initial data, we impose the following 

u Q G ^ r 2 D 2 (0\r c ; R 3 ) , [n ] h on S c , (3.7a) 
i. ei 2 (S];I 3 ) if > , (3.7b) 
z G L°°(T C ), < z a < 1 a.e. onr c , (3.7c) 
9 G L u (n) , 6» > a.e. in O, (3.7d) 

where w is as in (14. lbb below. 

Weak formulation. The energetic formulation associated with system ( 13.3b hinges on the following energy 
functional $, which is in fact the mechanical part of the free energy ( 12.13b . and on the dissipation metric 

31: 

$(u,z):=i / Ce(u):e(u) + HVe(u):Ve(w)d.T + / A '(H)+ / za ([u] ) + 7 [0 ,i] (z) dS, 
2 Mr c Jr c 

with the abbreviation a ([u] ) := ^A[u] • [m] - a ([u] ), (3.8) 



/ ai([u])|z|dS ifi < Oa.e. in T c 



R(u,z) := { J Fc VL J/l (3.9) 
-oo otherwise. 

We are now in the position of introducing the notion of weak solution to system ( 13.3b which shall be 
analysed throughout this paper. 

Definition 3.2 (Energetic solution of the adhesive contact problem). Given a quadruple of initial data 

(ito, Uo, zo, 8q) satisfying suitable conditions (cf. (13.7b later on), we call a triple (u, an energetic 
solution to the Cauchy problem for (the enthalpy reformulation of) system ( 13.3b if 

u G W h2 (0,T;W^(n\T c - : M. 3 )), (3.10a) 

u G W^°°(0,T;L 2 (ft;IR 3 )) if g > 0, (3.10b) 

2 G L°°(S C ) n SV([0, T];i x (r c )) , nonincreasing on [0,T]fora.a. x G T c , (3.10c) 

i?Gi r (0,T;W 1 ' r (O\r c )) n L^O.TjL 1 ^)) n S w »([0,T];M(fi)) 
?9 G W([0,T]; W 1 ' r '(0\r c )*) 



forany 1 < r < §, (3.10d) 



with r' denoting the conjugate exponent ^rj of r, and the triple (u, z, complies with: 
(i) (weak formulation of the) momentum inclusion, i.e.: 

[it] h onS c , and (3.11) 
q\ u{T) ■ (v(T)-u(T)) dx + I (De(u) + Ce(u)-MQ($)):e(v-u)- qu-(v-u) 



+ (lVe(u)+GVe(«)) Ne{v-u) dxdt + / a' ( [u] ) • {v-uj dSdt 

>q I u o -(v(0)-u(0))dx+ I F-(v-u)dxdt + [ f-(v-u)dSdt (3.12) 



for all win L 2 (0, T; W^{n\T c ; R 3 )) with [«] ^ on E c and, if q > 0, also in W ' (0, T; L 2 (0; R 3 )), 
(ii) total energy balance 



TKn(u(T)) + *(u(T),z(T)) + / 0(T)(dx) = T kin (u ) + $(u ,z ) + # dx 

Jq Jn 

+ / F-(v-u)dxdt + / f-(v-u)dSdt+ / ffdSdt, (3.13) 
Jq Js n is 

(iii) semistability for a. a. t G (0, T) 

Vz G i°°(r c ) : $(u(t), z(i)) < $(u(i), z) + %(u(t), z - z{t)), (3.14) 
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(iv) (weak formulation of the) enthalpy equation: 

w(T)0(T)(6x) + [ X(e(u),i})V$-Vw-tfwdxdt+ [ r]({u] , z) {&(■&)] {wj dSdt 
Jq Jt, c 

(Oe(u):e(u) - 0(i?)B:e(u) + GVe(u):Ve(u)) w dxdt 

-ai([u})z(dSdt) + I gwdSdt+ [ d w(0)dx (3.15) 



Hr c +Hr c 



for all w € C([0,T};W iy (n\r c ))nW 1 ' r ' {0,T; L r ' (ty), whereto := h (6 ), and i?(T) and z are 
considered as measures on and S c , respectively, 
(v) and the remaining initial conditions (beside u(0) = tto, already involved in (13.121 1). i.e. 

it(0) = ito a.e. in f2, z(0) = zq a.e. inr c , i?(0) = t?o a.e. in f2. (3.16) 

Remark 3.3 (The weak formulation ( 13.121 ) of the momentum inclusion). In order to (partially) justify 
(13. 12b and its link with the classical formulation ( 12.6b .c.d.f-k) of the (boundary-value problem for the) 
momentum inclusion, we may observe that, upon multiplying ( I2.6ab by v— u (with v an admissible test 
function in the sense of Definition l3.2l i and integrating on Q, one has to deal with the term 

— / div a ■ (v— u) dxdt = — / / diver KV -(i>— u) dxdt + / / drv 2 \)-(v— u) dxdt, 
jq Jo Jn Jo Jn 

where <r KV is a placeholder for the "Kelvin- Voigt" stress H)e(u)+Ce(u)— M8. The treatment of the first 
integral term on the right-hand side involves a standard integration by parts. As for the second one, let us 
observe (neglecting time-integration and integrating by parts twice, with the zero Dirichlet condition on 
T D ), that 



div 2 t) ■ (v— u) dx = — j (divh):V(w— u) dx — / (divfj):V(u— u) dx 
n Jn+ Jn- 



(divt)):((v-u)(E)is)dS + / (divh) + :((i>+-u + )(g)^c)dS 
(divh) _ :((v _ -M _ ) ® v c ) dS 

t)':V 2 (v-u)dx+ / t)':V 2 (v~u)dx 

Jn- 

+ I(T n ,u, v, v) + I(T c ,u + , v + ,v c ) - I(T C , u^^v^^c), (3.17) 

where we have used the short-hand notation I(T, u, v, v) := Jj,(divh):((5— u) <£> v)— h:(V(v— u) <E) v) dS. 
Then, the calculations developed in l32l 2nd ed., Sect. 2. 4.4], l36l and based on the decomposition V« = 



yield the following formula 

I(f,u,v,v) = [ ((divh)-!>+div s (h^)-(div s z>)((]:(£>(X)z>))) -(v-u) - " dS, 

if ov 

which we plug in ( 13.17b . Then, we combine the resulting integrals on T N and T c with the integrals derived 
from the by-part integration of J„ diva KV -(v—u) dx, rely on the boundary conditions d2.6H ,f-k). take into 

account the enthalpy transformation (13. Il l, and finally use that div 2 f) ■ (v—u) = f):Ve(u— u), since t) is 
symmetric, being so G and H. In this way, we obtain the second and the third term on the left-hand side of 
(13 .12b - The remaining terms either follow from an integration by parts in time, or are trivial. 

Remark 3.4 (The "weak" formulation of the flow rule ( 12.10b ). In J22][24]|26), a global stability condi- 
tion combined with energy conservation was shown to provide the correct "weak" formulation of rate- 
independent flow rules [22 24-26 1. Here, the concept of energy-preserving solutions (i.e. of energetic 
solutions) is crucial for mathematically treating the full thermodynamics, cf. Step 5 in Section|6]below. We 
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point out that d3.13l l is the integrated version of the total energy balance (12.20b . Here the energy conser- 
vation involves also the mechanical equilibrium ( 13.121 i. and the semistability ( 13.14b plays the role of the 
global stability condition of [22 24 26). We refer to l35l Prop. 3.2] for some justification of the energetic- 
solution concept in the framework of general thermomechanical rate-independent processes. In general, 
energetic solutions may exhibit unphysical jumps, but this does not occur if the driving energy ■) is 
convex, as it is indeed the case (13 - 8b considered here. Then there is also a close link to the conventional 
weak definition of the flow rule ( 12.10b . see 



Remark 3.5 (The weak formulation (13.15b of the enthalpy equation). A few comments on the first term 
on the left-hand side and on the second term on the right-hand side of ( 13.151 ) are in order. First, since 
i? G BV([0,T\; W iy (n\r c )*), then for all t G [0,T] one has well-defined as an element of 
W 1 ' r '(n\r c )*. Combining this with the fact that i? G L°°(Q,T; L 1 ^)) one sees that even G M(H) 
for all t G [0, T]. However, note that the function t M> -d(t) may jump. Second, let us observe that due 
to (13 .1 Oct . z is a negative Radon measure on Sc. Since we shall impose that the function oi : M. 3 — > M 
is continuous (cf. d4.1hb ). and since the map (t, x) H> lu(t, cc)J is also continuous because of (13. lOab and 
( 13.4b . it turns out that (t, x) i-> ai(fu(t, x)J) is a continuous function. Thus, ai([it])i is a well-defined 
measure on E c . 

Remark 3.6 (Mechanical energy equality). Subtracting ( 13.151 ) tested by 1 from ( 13.131 ) reveals that energetic 
solutions comply with the mechanical energy equality: 



T kin(«( T )) + ®( U ( T ), Z ( T )) + J ^Be(u):e(w)+GVe(w):Ve(w)J dxdt + Varj>(w, z\ [0, T]) 

= T kin("o) + ${u ,z ) + f F-u + S(i})M:e(u)dxdt + [ f-udSdt, (3.18) 

Jo JS» 



where we have used the notation 

Var^u, *;[<!,*,,]) := [' [_ oi([u])|«|(dSdt) for [t u t 2 ] C [0,T]. (3.19) 

4. Main result 

We now enlist our conditions on the functions c v , K, r], do, a\, and the loading. 
Assumptions. We suppose that 

c v : [0, +oo) — > (0, +oo) continuous, (4.1a) 

3wi > u > f, ci > c > 0V6» G (0,+oo) : c^l+Oy- 1 < c v (6) < , (4.1b) 

K : ]R 3x3 xR-> E 3x3 is bounded, continuous, and (4.1c) 

inf X(e,tf)£:£ =: k > 0. (4. Id) 

(e,i9,{)eRs y x „?xRxR 3 , |{|=1 

We also require that i](x, v, ■) is a non-negative affine function of the delamination parameter z G [0, 1], 
and, following OTI . we assume that 

rj(x, v, z) = rji(x, v)z + r] (x, v) for 771,770 : r c xR 3 — > [0, +00) Caratheodory s.t. 

3C V >0 V(x,v)eT c xR 3 : \r] (x, v)\ + \ m {x, v)\ < C v {\v\ 4 / 3 + 1); 

As for the functions gq : R 3 — > R and a\ : R 3 — > R, we suppose that 

aoeCyjR), 3C a „,C: o >0V U £K 3 : |a (v)| < C ao H + C' ao (4.1f) 

the map u h->- -Au-u — ao(u) =: ao(u) is convex, (4.1g) 

djGCft 3 ;!), 3 O ai > Vij G R 3 : oi(v) > C ai > 0. (4.1h) 
Remark 4.1. Let us comment on conditions ( |4.1b . First of all, it is immediate to deduce from ( |4.1bb that 

3 C e \ C 2 e > Vttj G [0, +00) : C e V 1/aJl - 1) < B(iu) < C^ttj 1 /" + 1) . (4.2) 



(4.1e) 
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It obviously follows from (14. let that 

3C X >0 V£, (eR 3 : |3C(e,i?^:C|<C x ^||C|. (4.3) 

Moreover, let us observe that the functional u i-> $(it, ar) is convex thanks to ( |4.1g[ i. Note that ao itself need 
not be concave, and the possible violation of concavity depends on the positive-definiteness of A. Actually, 
we could even allow for bigger violation (namely for ao semi-concave), if the discretization scheme were 
slightly modified, like for example in 1321 2nd ed., Rem. 8.2.4]. However, we have chosen not to explore 
this option, since in real-world applications A is large. 

Theorem 4.2 (Existence for the adhesive contact problem). Let us assume (13. 6b , (14. U . ( 13.7b and 

(i) if g = 0, suppose also 

FeW u (0,r;L 6 / 5 (S];i 3 )), 

jf 2 (dn + n r D ) > o, jf 2 {dn_ n r D ) > o, 

where Jt? 2 denotes the two-dimensional Hausdorff measure, or 

(ii) if g > 0, suppose also that 

K{x) is a linear subspace of R 3 for a.a. x <G T c . (4.5) 

Then, there exists an energetic solution (u, z, to the adhesive contact problem (in the sense of Definition 
13.2b . with the additional regularity 

u€W h2 (0,T;W^ 2 (n\r c ;R 3 )*) if g > 0. (4.6) 

Furthermore, in both cases g > and g = 0, the positivity of the initial temperature 

inf 6» =: 0* > (4.7) 

implies inf^^gg 6 = m£n x \ e Q 0(i?(i, x)) > 0; in particular, 8 is a.e. positive on Q. 

Remark 4.3. The analytical reason why in the presence of inertial terms in the momentum equation we 
need to restrict to "linear" contact conditions on T c is ultimately that, if g > 0, only ( 14.51 ) makes it possible 
to test the (weak formulation of the) momentum equation by the velocity ii. This is needed for obtaining 
the mechanical energy equality ( 13.181 l, which in turn is a crucial step in the proof of Theorem |4.2| 

In what follows, we shall denote by the symbols C, C most of the (positive) constants occurring in calcu- 
lations and estimates. 



(4.4a) 
(4.4b) 



5. Semi-implicit time discretization 

We perform a semi-implicit time-discretization using an equidistant partition of [0, T], with time-step 
t > and nodes := kr, k = 0, . . . , K T . Hereafter, given any sequence {(/)•' }j>i, we will denote the 
backward difference operator and its iteration by, respectively, 

„ t=£L\ D y :. D,(D,^) J±^±fl ' (5 „ 

T T 

We approximate the data F, f by local means, i.e. setting for all k = 1, . . . , K T 

F k ; := - t F(s) ds , f k T := - t f{s) ds . 

Furthermore, we approximate g by suitably constructed discrete data {g^}^!-, C H 1 ^ 2 (dfl)* such that 
( 15.101 1 below holds, and the initial datum m by a sequence {w , T } C Wp'J (n\T c ; R 3 ) (with 7 > max{4, ^5-} 
as assumed in Problem [5. U such that 

Urn ^/r||Ve(it , T )||L-v(n;R3x3x3) = 0, u , T -> u Q in W 2 ' 2 (n; R 3 ) as r ->• 0. (5.2) 

We are now in the position of formulating the time-discrete problem, which we again write in the 
classical formulation for notational simplicity. 
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Problem 5.1. Let 7 > max{4, ;jtt}. Given 



Mo, 



u T 1 = it , T - TM , z° = z Q , d° r =-d , 



(5.3) 



find {(u k ,^ k ,z k )}^Z 1 fulfilling, for k = l,...,K T , the recursive scheme consisting of the discrete mo- 
mentum equation in ri\r c : 



gB 2 t u k 



^Be(Dj^) + Ce(u*)-B0(0*) +r|e(w^)| 7 2 e{u k T ) - divh^ = F k 



with 



+ T|Ve(^)| 7 2 I) Ve(u£)+GV(e(D t u*)) 



in ^\r c 



where I : 



^3x3x3 ^3x3x3 d eno t es m e 6th-order identity tensor, with the boundary conditions 

on r D , 



(De(D t u k ) + Ce(u k ) - 0(^)1 



+ r|e(^)| 7 ^(O-divh^)^ 
b*:(i/®J/) = 
and the conditions on the contact boundary 



on r N , 
on r D u r N , 



pJe(D t u*) + Ce(u*)-0(0*)B + T|e(i4)| 7 2 e(it^) - div(h*)K 

-div 8 ([h*]i/ c ) = 0, 

(h*) + :(ifcg>iO = {t) k ;)-:(v c ®v c ) = 0, 

zJ-Vo ([«*]) + dI K (H}) + (Be(F> t u k ) + Ce(u*)-0(0*)B 

+T\e{u k )\~ < ~ 2 e(u k ) - divf^W - div s (h* • u c ) = 0, 



> on r r 



further, the discrete enthalpy equation: 

D t l?*-div(3C(^ J e(«*))V^) =De(D t «*):e(Dt«*) 

- 0(i?*)B:e(D t tt*) + GVe(D t tt*):Ve(D t u*) 
with the boundary conditions 

(X{d k ,e{u k ))Vd k )- V = g k onr D ur N , 
and the conditions on the contact boundary 

i(3C(^,e(^))V^|+ + 3C(^,e(^))V^|F c )^ c 

+V(H- 1 } ,z k ) [0(^)] =0, 
iX(d k ,e(u k ))\7d k p c = -oi([«*])D t «* 

and also the discrete flow rule for the delamination parameter 

9„J(z r fe - 1 ;z T fe ) + a ([^])-a 1 ([^]) 3 on T c , 
where ST^* -1 ; •) : M — > [0, +00] is the convex functional 



in n\r c 



on rv 



y(^- 1 ;«) = J ( 



(-00,0] I 



J[o,i](w)- 



(5.4a) 
(5.4b) 

(5.4c) 
(5.4d) 



(5.4e) 



(5.4f) 



(5.4g) 



(5.4h) 



(5.4i) 



(5.5) 



In the last condition in ( |5.4et , traces of the overall stress either from Sl + or from Q_ can be considered 
with the same effect, thanks to the first boundary condition in d5.4et . 
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Remark 5.2. Let us highlight the main features of the time-discrete scheme (15.41 ). 

First, the discrete version ( I5.4f1 i of the enthalpy equation is fully implicit, and in particular on the right- 
hand side the term 0(^)B:e(DtU^J appears, instead of Q(d k \~ 1 )M:e{ptu k \) . This is crucial to obtain the 
positivity of the temperature, i.e. & k > a.e. in il, cf. Lemma [5~4l Notice that the term B0(i9^ _1 ) occurs 
on the left-hand side of ( I5.4at : therefore, ( 15.4al i and (I5.4f1 i are coupled. 

Second, the boundary conditions ( I5.4hb on T c for the discrete enthalpy equation involve z k . In this 
way, (15 .4fb is coupled with ( I5.4il i. hence the whole system is coupled. Nonetheless, the mechanical part of 
system ( 15.4b (viz., d5.4b -e.i)) can be reformulated in terms of the subdifferential inclusion 

5T kin (D^) + d ( . ti) S(u k T ; e(B t u k T ),D t z k ) + d M *(u k , 0*) 3 L k T (5.6) 

(cf. ( 12.16b ). i.e. it is semi-implicit w.r.t. the variable z. This will allow for some simplifications in the a 
priori estimates, see Lemma lBTol 

Third, we have added the term rdiv 2 (T|Ve(it^)| 7 " 2 Ve(i4)) - Tdiv(T\e(u^' 2 e(u k )) to the mo- 
mentum equation in the bulk and to the corresponding boundary/contact conditions, too. Its role is to 
compensate the quadratic growth of the right-hand side of the enthalpy equation (15.4fb when 7 is chosen 
large enough, cf. the proof of Lemma 15~4l Being premultiplied by the factor r, this higher-homogeneity 
regularization will vanish when passing r — > 0. Because of this term, we also need to regularize the initial 
condition uq in ( 15.31 , cf. (15.2b . 

Remark 5.3. The time-discrete scheme ( 15.4b is simpler than the one devised in PP . where the first term 
on the right-hand side of (I5.4fb was multiplied by the coefficient (1— \fr), and additional terms (featuring 
monotone functions of z and \u\ ) were added to the discrete flow rule and to the boundary conditions on 
r c for u. Such terms were used in the derivation of the discrete a priori estimates (in particular, of the 
discrete energy inequalities) via auxiliary minimization problems. Instead, here we adopt a more direct 
approach in the proof of the discrete mechanical and total energy inequalities, cf. ( 15.16b and ( 15.13b below. 
Indeed, we strongly rely on the semi-implicit character of ( 15.61 ) and on the convexity of o.q. 

Lemma 5.4 (Existence of weak solutions to Problem l5.1b . Under the assumptions of Theorem 14.21 for 
every k = 1,...,K T , there exists a triple (u k ,z k ,$ k ) £ W^{fl\T c ; M 3 ) x L°°(r c ) x W^^Ic), 
fulfilling the weak formulation of the boundary value problem (I5.4ab - (l5.4il ). Moreover, d k > a.e. in fi. 
If, in addition, (14. 7b holds, then there exists some constant x* > such that, for sufficiently small t, 

i?r > X* > a.e. in d for every k = 1, K T . (5.7) 

Sketch of the proof. We may argue along the very same lines as in the proof of OTl Lemma 7.4]. Indeed, 
the existence of a weak solution to Problem 15.11 follows from the theory of pseudomonotone set- valued 
operators (see e.g. 11321 Chap. 2]), and in particular from Leray-Lions type theorems, like ll32l Chap. 5, 
Cor. 5.17]. To apply such results, one has to verify the strict monotonicity of the main part of the ellip- 
tic operator, involved in the weak formulation of problem ( 15.4ab -( f5~4il i. One has also to show that this 
operator is coercive w.r.t. the norm of W^(n\T c ] R 3 ) x L°°(r c ) x W 1 ' 2 (ri\T c ). For this, the term 

Tdiv 2 (|Ve(M^)| 7-2 Ve(i4)) — rdiv(|e(u£)| 7_2 e(u£)) on the left-hand side of d5.4ab plays a crucial role, 
in that it counteracts the quadratic nonlinearities in e(u k ) and in Ve(u£) on the right-hand side of ( 15.4ft : 
for this , we need 7 > max{4, -j^r }. All the calculations for proving this strict monotonicity and coerciv- 
ity in the present setting, and also for obtaining the strict positivity ( 14.7b . are very similar to those carried 
out in the proof of Oil Lemma 7.4]. Therefore, we prefer to omit all details and directly refer the reader 

to on. □ 

Remark 5.5. Note that the actual discrete version of the flow rule ( 12.10b for the delamination parameter 
would be 

5J ( _ 00>0 ](D t ^)+a ([«*])-oi([«*]) + W [ o,i](^)9 onr c , (5.8) 
which in fact yields a solution to d5.4ib , since dl^^ m (p t z k ) + 9/fo ) i](^) C d v 3 r (z k ~ 1 ; z k ), while the 
converse inclusion in general does not hold. The reason why we have replaced ( 15.8b by ( 15.4ib is that, 
differently from d5.8l ). the differential inclusion d5.4ib features only one nonsmooth unbounded operator. 
Hence, we are entitled to directly apply the aforementioned l32l Chap. 5, Cor. 5.17] to prove existence 
of solutions to Problem [5. II In turn, in view of the analysis which shall be developed later on (cf. Lemma 
15.61 ), it is actually sufficient to solve d5.4ib in place of d5.8b . 
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Approximate solutions. For r > fixed, the left-continuous and right-continuous piecewise constant, 
and the piecewise linear interpolants of the discrete solutions {v%}?I-, are respectively the functions Hj- : 
(0,T) -> W^ D 7 (ft\r c ;IR 3 ), u T ■ (0,T) -> W^(n\T c ; R 3 ), and : (0,T) -> Wp^ 7 ( il\ T c ; R 3 ) defined 
by = u*, = u* _1 , Ur(i) = — u k + ^^u^ -1 fort e (t* -1 ,i*i]. In the same way, we shall 
denote by i9 T , t? T , and z T , the piecewise constant interpolants of the elements {'& k }^Z 1 an( ^ { z r}£i' an d 
the related piecewise linear interpolants by i? T and z T . Furthermore, we shall use the notation t T and t T for 
the left-continuous and right-continuous piecewise constant interpolants associated with the partition, i.e. 
lr(t) = t k if t*' 1 <t<t k and| T (» = if** -1 <t<t k . 

We shall also consider the interpolants F r , / T , and f T , of the A' T -tuples {-F^ {fr}k=i- ^ n y i ew 
of ( I3.6ab — ( f3.6bb and ( I4.4al ). the following estimates and strong convergences hold as r — s- 0: 

= f jnL 1 (0,T;L 2 (n;R 3 )) if e >0, 

r \ m 17(0, T;L 6 / 5 (f2;R 3 )) for all l<p<oc if p = 0; V ; 

3C > Vr > : ||/ T ||L~(UT;L 4 / 3 (r N ;R 3 )) < ^I!/IIl~(o,T;L 4 / 3 (T n ;R 3 )) > 
7 r -> / inL p (0,T;£ 4/3 (r N ;R 3 )) for all 1 <p<ooasr^0, 

3C > Vr > : ||/T||L!(UT;L 4 / 3 (r N ;R 3 )) < 2||/|| L i( ,T ; L 4 / 3 (r N; R 3 )) , 

3C > Vr>0 : ||-fV||Li(o,T;L 6 / 5 (n;R 3 )) ^ 2 ll- F llL 1 (0,T;L 6 / 5 (n ; ]R 3 )) i 



(5.9b) 



Finally, we shall construct the discrete data {g k }^l x C H 1 ^ 2 (dn)* in such a way that the related piecewise 
constant interpolants g T fulfill 

g T ^g ini x (E) as r ->0. (5.10) 

Using the interpolants so far introduced, we now state the discrete versions of the weak formulation (13.12l i 
of the momentum inclusion, the total energy balance ( 13.131 1, the semistability (13.141 1, the weak formulation 
(13.151 of the enthalpy equation. For the momentum inclusion, we introduce "discrete test functions", viz. 
AVtuples {v*}^ C W^ D 2 (fi\r c ;R 3 ), fulfilling \v k \ ^ on T c , and we denote by v T and v T their 
interpolants. Furthermore, referring to the definition ( 13.81 ) of $, we shall use the notation 

® T {u,z) := $(u,z) + - / |e(u)| 7 + |Ve(u)| 7 da;. (5.11) 
7 Jn\r c 

Hence, the approximate solutions (jL T ,u T ,d T ,^ T ,~z T ,v n -,'d T ,z T ) fulfill the discrete (weak) momentum 
inclusion 

(De(«r) + Ce(ur) - B6(tf T ) + -r|e(l^) | 7 '~ 2 e(TZ r )) :e{v T -Ur) 
+ ((H + r|Ve(u r )| 7 ~ 2 I)Ve(u r ) + GVe(i T )):Ve(iJ r -u r ) ) dxdt 



Z^dur] )-[v T -UrJ dxdt - / / 

Ji Jo 

+ / ^(T)-K(T)-u r (T))d.T 
Jo 



gii T (- — t)-(v t — iir) dxdt 



> / QU Q -(v r (T)-v T (T))dx + / F r -(?J r - %-) dxdt + / / T -(u r - %■) dSdt; (5.12) 

JO JQ JS N 

the discrete total energy inequality 

T]dn(ur(*)) +*r(Ur(*),3r(*)) + / r (t)dB < T kin («„) + $ T (u . T , z ) 

JO 

+ f $ dx+ f(! F T -ii T dx + [ J T -{wdS+ [ g T ds)ds; (5.13) 
Jo Jo \Jo Jr N Jao / 
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(see ( 12.15b for the definition of T kin ), the discrete semistability for a.a. t e (0, T) 

* r fc(t),2j.(t)) <$ T («r(t),«)+3l([^(t)],5-^.(t)) forallzei 00 ^); (5.14) 
the discrete (weak) enthalpy equation 

d T {T)w{T)dx + I X(e{u r ),d T )V^ T -Vw-Awdxdt+ /?y([w T ] ,z r )[e(^ T )] [wj dSdt 



7>e(ur):e(ur) - 6(t9 T )B:e(?V) + GVe(?V):Ve(?l r )J wdzdt 

- / a^lvir})^ c c dSdt+ / i!) o w{0)dx + / g T wdSdt. (5.15) 

with w qualified as in ( 13.15b . Inequality ( 15.12b can be obtained from d5.4b -e). by using a suitable discrete 
"by-part" summation formula, cf. (35] Formula (4.49)]. d5.4f -h). We now prove ( 15.13b and ( 15.141 ). 

Lemma 5.6 (Approximate energetics). Let g > 0. Under the assumptions of Theorem l4.21 for all r > the 
approximate solutions (%■ , u T , i? T , ~z T , Ur , $ r , Zr ) fulfill the following discrete mechanical energy inequality 

rt T (t)/ r f \ 

)e(ur):e(ur) +GVe(ii rr ):Ve(u r ) dx + J &([%■] , z T ) dSJds 

<T iin (u ) + ^ T (u , T ,z ) + f ( f e(^.)B:e(«r) da + / F r -u T dx + f J r -^ds]ds, 

Jo \Jn Jo Jr N / 

(5.16) 

(see ( I2.12bl ) for the definition of the discrete total energy inequality ( 15.13b . and the discrete semistabil- 
ity (EH. 

Proof. Preliminarily, we observe that ( I5.4ib is the Euler-Lagrange equation for the minimum problem 

z k G Argmin 2eL oo (r } / za ([«*])+ I?^ 1 ; z) dS, (5.17) 

where 5F(z* _1 ; z) is as in (15.5b . Therefore, we have 

< / z^([^])-^V(K])dS. (5.18) 



Jr c 

To prove ( |5.16b , we test the boundary-value problem J5.4a| >— d5.4et by u k —u k ~ l . We add to the resulting 
inequality the previously observed ( 15.18b . thus obtaining 

q [ B^u k -D t u k dx + [ De(Dt«*):e(D t «J)+GVe(D t u*):Ve(D t «*)da! 
Jn Jn 

+ f Cc(uJ):e(D t u*)+HVe(«*):Ve(D t «J)da: 
Jn 

+ t f |e(^)r- 2 e(^):e(D t ^)+|Ve(^)|T- 2 Ve(^):Ve(D t ^) dx 

+ / z^a'oilu^iD^j ~ (D t z k T )ao(lu k T j)dS - [ oi([«*] XD^" 1 ) dS 
=: /i + I a + J 3 + Zi + ^5 + h 



< T 



[ B0(^):e(D t u*) dx + r / F k -D t u k dx + t [ f k -D t u k dS. (5.19) 
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Now, we estimate the terms Ii, i = 1, . . . , 8. First of all, we observe that 

h>-!r l |D t ^| 2 dx- ± [ |D t ^- 1 | 2 dx=|D t / |D t ^| 2 cb. (5.20) 
2r J Q 2t J q 2 J n 

Clearly, upon summation I 2 yields the third summand on the right-hand side of (15.16b . whereas we observe 
that 



h > D t /_i(Ce(«*):e(«*)+HV(e(t^)):V(e(u*)))da:, 



n 

T , 



h > D t / -|e(u*)p r + -|Ve(u*)Pda; J 
in 7 7 



(5.21) 



where we have used elementary convex-analysis inequalities. Now, it follows from the convexity of ao, cf. 
( |4.1g) , that 

/ zk-^Wflyptv*] dS>- ( z^a (lu k T ])dS-- ( ^([^DdS. 

Hence, taking into account the cancellation of the term J r z k .~ 1 ao(lu k \j), for Is we conclude the following 
inequality 

7 5 >D t / 2*a ([«*])dS'. (5.22) 
Jr c 

Combining (I5.19b - d5.221 i, rearranging terms and multiplying by r, we obtain 

/ ||D f ^| 2 + rDe(D t ^):e(D t ^) + rGVe(D t uJ):Ve(D t uJ) + -Ce(u*):e(u*) 
Jn 2 2 

+ -| e (^)r + -|Ve(^)r + iMV(e(^));V(e(^))d a; + f z r fc a ([u*])dS 
7 7 A Jr c 

< I e \V t u k T - 1 \ 2 + \ I Ce «- 1 ):e(^- 1 ) + -| e (^- 1 )P + I|Ve(^- 1 )r 
Jn 2 2 Jn 7 7 

+ ±UV(e(u k ;- 1 )y.V(e(u k ;- 1 ))dx+ f z^aoilu^dS 

+ t I Ke(^):e(D t u*)da; + T / F k -T> t u h ; dx + t [ j k -D t u k T dS. (5.23) 
Jn Jn Jr N 

Summing over the index k, we conclude ( 15.161 ). 

In order to obtain ( 15.13b for a fixed t <E (0, T), we test (15.15b . integrated on the time-interval (0, tr(t)), 
by 1, and add the resulting relation to the mechanical energy equality ( 15.16b . Taking into account all 
cancellations, we immediately conclude ( 15.13b . 

Eventually, from (15.18b and the degree-1 homogeneity of ^([uj;]] , •), it also follows that 

$ r {u k T ,z k T )<$ r {u k T ,~z)- f a 1 ([^])(5-z^ 1 )d5- f ^([uJilXz*- 1 -*^ 

Jr c Jr c 

<$ T (u k ,z)- f oi( [«*])(«-«*) d5 = * T (u*, 2) +3l([uj!], 2-^) 

JTr 



(5.24) 



for all z < z k a.e. on T c . This is the discrete version of (15.14b . □ 
We conclude this section with a result collecting all the a priori estimates on the approximate solutions. 
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Lemma 5.7 (A priori estimates). Under the assumptions of Theorem 14.21 for all g > and t > 0, the 
approximate solutions (ur, <& T ,z T ,v T ,'d T ,z r ) satisfy 



\\ Ur L°=(0,T;Wp D 2 (O;R 3 )) — ^° ' 


(5.25a) 


K" llw 1 - 2 (0,T;W r ^ 2 (n;R 3 )) - ^° ' 


(5.25b) 


1/211 II < - a 
\\^\\w 1 ' oa (0,T;L*(n-B. 3 )) — ' 


(5.25c) 


M_ || < So 

ll Ur llL~(0,T;Vy 2 D "'(a;R 3 )) - "^9' 


(5.25d) 


\\ Jr IIl°°(e c ) - S ° ' 


(5.25e) 


\\ Zt BVdO.T];! 1 ^)) - 


(5.25f) 


lr T lli oo (0,T;i I (n)) — 


(5.25g) 


W4 L r(0,T;W^ m ^ S r for ^ 1 < ^ < |, 


(5.25h) 


r T IL 1 (o,T ; w rl .'-'(n)*) - 


(5.25i) 


8 sy([o,T] ; w ] 2 D ' T (a;R 3 )*) — ^ ' 


(5.25j) 



for some constants So > and S r > independent of r. Estimates ( 15.25b . f) hold for 2v as well, and so 
do estimates ( 15.25fc .h) for d T . 

Proof. We only sketch the calculations for proving ( 15.25b . since the argument closely follows the proof 
of OTl Lemma 7.7], to which we shall systematically refer. 

First of all, we use the "discrete total energy" balance ( 15.13b . Clearly, the first summand on the left-hand 
side provides a bound for the quantity g 1//2 ||itr|| w ,i iOO ^ T-L 2 (n-n?))' Secondly, we observe that (cf. (13. 8k 

•mux 

$ r (u, z)>c( K J \e{u)\ 2 + |Ve(u)| 2 + r\e(u)p + r|Ve(u)| 7 dx 




where we have used the positive-definiteness of A, C, and H, and the growth condition ( 14. 1 ft . to derive 
the first inequality. The second estimate ensues from Korn's inequality, and from absorbing the term 
J r z\fuj | da; into J r ^ z\fuj | 2 dx, since z 6 [0, 1] a.e. on T c . Therefore, the second term on the left-hand 
side of ( 15.13b estimates \\ur(t)\\ ^ (f2 . R3) and T||^(i)||^/2, 7(a . R 3) uniformly w.r.t. t G [0,T], Thirdly, 

■d T > a.e. in O thanks to Lemma [5~4l hence the third term estimates \\fl T \ \ (o,T:L~>- (n)) ■ To deal with the 
right-hand side of ( 15.131 ). we use ( 13.7b . ( 15.21 ) and, for the last integral term, ( 15.91 ) and ( 15.101 ). arguing in the 
very same way as in the proof of [31, Lemma 7.7]. We conclude applying the discrete Gronwall lemma, 
and thus obtain estimates ( 15.25al ). ( 15.25cb . and ( |5.25g[ ). Since z r e [0, 1] a.e. on E c , we obviously have 

Secondly, again arguing as for EH Prop. 4.2] and EH Lemma 7.7], we make use of the technique by 
Boccardo and Gallouet |3i] with the simplification devised in JTJ], and we test the heat equation ( 15.151 ) by 
7r(# T ), where n : [0, +00) — > [0, 1] is the map w i-> n(w) = 1 — (1 } w) ^ , for <r > 0. Since n is Lipschitz 
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continuous, 7r($ T ) € W^ 1,2 (r2\r c ) is an admissible test function. We thus have 

|2 



k j |W T | — d d /" ^YgfaN # r )V^ T -V7r(^ T )dxdi 

77([Mr],^)[e(^r)][^ r )] dSdt+ I n(MT,-))dx 



£(0o)da: + C'(J|Be(^):e(^)|| L i ( Q ) + ||GVe(u r ):Ve(u r )|| L i (Q) 
+ \\Q@ T )M:e(ur)\\ LHQ) + ||Ci([^] , ^Olk^c)) + W&hw < 5 - 27 ) 



< 



where 7? is the primitive function of tt such that 7r(0) = 0. Note that inequality ( 15.27b follows from ( 14.1db . 
the fact that 77 (u T ,z T )[{ 6 ($,-)]] [[7r($ r )]] > a.e. in S c (by the positivity of t] and the monotonicity of 
and 7r), from the "discrete chain rule" l35l Formula (4.30)] for tt, and from the fact < 7r(^ T ) < 1 a.e. in 
ft. Combining d5.271 i and performing the very same calculations as in the proof of 051 Prop. 4.2], with the 
Gagliardo-Nirenberg inequality we find for all 1 < r < 5/4, that 

ll V ^ll^(Q;R3) < a(l + \\De(ur):e(ur)\\LHQ) + ||<GVe(^):Ve(?V)|Ui (Q) 

+ ||e(^)l:e(^)|| iI(Q) +||Ci([u r ],iv)|U 1 ( Sc )) (5.28) 

for some positive constant C r , depending on r and also on the function 77, cf. (14. let . 

Then, we multiply d5.281 l by a constant p\ > and add it to ( 15.161 ) (in which we set t = T). Now, by 
positive-definiteness of D and G, the third term on the left-hand side of ( 15.161 ) is bounded from below by 
c (ll e (^)ll|2(Q ; R3x3) + ||Ve(w T )|li2 (Q;R 3x3x3 ) ), and it controls HCife], 2t)||l 1 (s c )- Thus, we choose pi 
small enough in such a way as to absorb the second, the third, and the fifth term on the right-hand side 
of (15.28b into the left-hand side of (15.161 ). Hence, we find 

C (ll e (Wr)llL2(Q;R3x3 ) + ||Ve(u r )||i2(Q; H 3X3X3)) + Pi) || & ( [VJ , *T ) \\v (E C ) + Pi \ | V#r 1 1 \ r (Q;R 3 ) 

< T kin (u . T ) + <$> T (u , T , Z , T ) + / F r -{l rr dxdt+ J % -Vr dSdt + (p 5 C r +l) || 0(^.)Be(tir) 1 1 ^ {Q) . 

(5.29) 

The first two summands on the right-hand side of ( 15.291 ) are estimated in view of ( 13.71 ) and ( 15.21 ). Using 



( 15.9b , we handle the terms f Q T J" n F T -Ur dxdt and / Q T J r f T -Ur dSdt in the very same way as in the proof 



of OTl Lemma 7.7]. Finally, we use 

( Pl a+l)||e(^0B:e(^)||L 1( Q) <P2||e(^)||i2 (Q;R 3x3 )) +C' P2 ||e(^)l| 2 L2 (Q) 

< P2 l|e(iV)|||2 (Q . IR 3x 3)) + C P2 (p T ||^( Q )+ !) 

< P2 ||e(^)|li2 (Q;R 3x3 )) +^3^ ||V^ T ||^ r(0;K3) di + Cp 3 , (5.30) 

where the last inequality can be proved, via the Gagliardo-Nirenberg inequality, by developing the same 
calculations as throughout l35l Formulae (4.39)-(4.43)], and using the restriction on oj in ( 14. lbb and the 
previously proved bound for ||#r||i~(o,T;£ 1 (f2)V Then, we plug ( 15.30b into ( 15.291 ), and choose pi and 
P3 in such a way as to absorb the terms ||e(u,-)||i2(g. R 3x3) and || Vi? T ||^/Q. ffi3 N into the left-hand side 
of ( 15. 29b . Thus, we conclude estimate d5.25bb . We also get an estimate for HCidf^-IK ■Zt)IIl 1 (i: c )« which 
yields ( I5.25fb . since a\ is bounded from below, cf. ( I4.1hb . Furthermore, we also obtain a bound for Vt9 t 
in L r (Q;M. 3 ). Combining the latter information with the estimate for d T in L°°(0, T; L 1 (i7)), we in- 
fer 45223. 

Estimate (I5.25ib follows from a comparison in (15.15b . and the related calculations are a trivial adaptation 
of the ones in the proof of DTI Lemma 7.7]. 

Finally, for (5.25j 1 we use that iL is a measure on [0, T], supported at the jumps of Ur, and we estimate 
£RllM(0,T;W^(fi\r c ;K3)*) by comparison in (ETJ. □ 
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6. Proof of Theorem |4.2| 

In what follows, we develop the proof of the passage to the limit in the time-discrete scheme as r — > 
unifying the cases g > and g = 0; we shall take a sequence of time-steps, i.e. we understand (r) as 
countable family of indexes with the accumulation point 0. 

For the reader's convenience, we briefly describe the strategy. After a careful selection of converging 
subsequences is made in Step below, by passing to the limit as r — >• in ( 15.12b we will obtain the 
weak formulation ( 13.12b of the momentum inclusion. Then, we will proceed to proving the semistability 
condition ( 13.141 ). hence the total energy inequality by lower semicontinuity arguments. By the same tokens 
we will also obtain the mechanical energy inequality. We will then show that the latter in fact holds as 
an equality, by combining a chain rule-type argument (cf. ( 16.191 )). with a test of ( 13.121 ) by u. To perform 
the latter, it will be essential for u to have the regularity ( 16.28b . This motivates the dissipative contribution 
GVe(ii) to the hyperstress. Hence, we will exploit the mechanical energy equality to conclude, via a 
suitable comparison argument, the convergence of the quadratic terms in the right-hand side of ( 15.151 ). This 
will allow us to pass to the limit, and conclude the weak formulation ( 13.15b of the enthalpy equation and, 
ultimately, also the total energy balance ( 13. 131 ). 

Step 0: selection of convergent subsequences. First of all, it follows from estimates (I5.25bl ). ( I5.25cl ). and 
( |5.25j| i, from the Banach selection principle, the infinite-dimensional Ascoli and the Aubin-Lions theorems 
(see, e.g., [41 Thm. 5, Cor. 4]), that there exist a (not relabeled) sequence r —> and a limit function 
u G W 1 ' 2 ^, T; W^(n\T c ; R 3 )) such that the following convergences hold as r -> 0: 

Ur^u mW h2 (0,T;W^(n\T c ;R 3 )), (6.1a) 

Ur^u inC([0,T]; W^^ e ' 2 (fi\r c ;R 3 )) Ve e (0,2], (6.1b) 

gur^-gu inW 1,oo (0,T;L 2 (il',R 3 )). (6.1c) 

Estimate ( 5.25j| l and a generalization of the Aubin-Lions theorem to the case of time derivatives as measures 
(cf. e.g. J32] Cor. 7.9]) also yield that u e BV([0, T]; W^{n\T c ; R 3 )*) and that 

gitr^gu ini 2 (0,T;W^^ £ ' 2 (fi\r c ;R 3 ))forallee (0,2]. (6.1d) 

Moreover, a generalization of Helly's principle (see |2| as well as l26l Thm. 6.1]) implies that Ur(t)^ u(t) 
in Wp D 7 (f2\r c ; R 3 )* for all t £ [0, T]. In view of estimate ( I5.25cb . with an elementary compactness 
argument we conclude 

gUr(t)^Qu(t) ini 2 (ft;R 3 ) for alU e [0,T]. (6.1e) 

Next, we observe that 

||^r-Ur|| L ^ ( o, T ;W r 2 D 2 (o\r c ;B3)) < rV2 li «r || L 2 (0)T . w p* (n \r c ;R3)) < S t 1/2 ^0 as T -> 0. (6.1f) 

Therefore, estimate d5.25ab and d6.lab - d6.lbb yield for all e S (0, 1] 

Hr^u mL oo (0,T;H§j I (n\r c ;R 8 )) ) (6.1g) 

Ur^u mL°°(0,T;W 2 ^ 2 {n\T c ;R 3 )), (6.1h) 

Ur(t)->u(t) iniy 2 D ^ 2 (fi\r c ;R 3 ))foraine [0,T]. (6.1i) 

Taking into account the compact embedding (13.4b . from ( 16. lbb and ( 16. lhb we deduce respectively 

[ur] -> [«] inC(S c ;R 3 ), (6.1j) 

[ur] -> [«] inL°°(I] c ;R 3 ). (6.1k) 

In fact, the above convergences are also in C([0, T];C(T C ; R 3 )) and in L°°(0, T; L°°(r c ; R 3 )), respec- 
tively. Convergences (16.1^ -i.k) hold for u T , too. Also note that, in view of d5.25dl ). we have 

^|||e(^)l 7 " 2 e(^)|L,/„-i, (Q:K 3x3 ) <5 rVT^0 and 1 

r|||Ve(u r )r- 2 Ve(u r )|| iV( ,_ 1)(0;R3>< 3 X 3 ) < S r^^ ] as r ^ U ' ^ 

Estimates d5.25e| ) and d5.25fT ). and the very same compactness arguments as in OTl Sec. 8] (based 
on |26l Thm. 6.1, Prop. 6.2]), also guarantee that there exists a function z e L°°(T, C ) H BV([Q,T}; Z) 
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(where Z is any reflexive space such that L 1 (r c ) C Z with a continuous embedding), such that, possibly 
along a subsequence, 

z T z mL°°(T, c ), z T (t) z(t) in L°°(r c ) for all t G [0, T\. (6.2) 
We now prove that 



Var 3? (w, z; [s, t]) < lirninf / / & , z T ) dSdr for all < s < t < T. (6.3) 



Indeed, 



( 1 (lu r ],i T )dSdr= / / a^lurj^ZrldSdr 



ai ({u T ])\z T \dSdr+ / / ( ai ([u r ])-a 1 ([« T ]))|i r |d5dr. (6.4) 
/r c Jr c 

Now, from (I6.1| ,k) and the continuity of a\ it follows 

at({ur]) - ai([ur]) -> Q mL°°(0,T;C(T c )). (6.5) 

Since (z T ) T >o is bounded in L 1 (S c ), we then conclude that the second term on the right-hand side of ( 16.4b 
tends to zero as r — > 0. To pass to the limit in the first term, we use (6.1j I and again the continuity of a\. 
Since \z T \ — > \z\ weakly* in the sense of measures on S c , we conclude that 

ai([u T ] )|ir| — > ai(M )|i| weakly* in the sense of measures on S c , (6.6) 

where \z\ denotes the variation of the measure 'z\ in fact, here simply \z\ = —z, since 'z < 0. Then, ( 16.3b 
follows. Taking into account that a\ is bounded from below by (14. lhb . (16.3b and the definition ( 13.191 ) of 
Var* imply z G BV([0,T}; L\T C )). 

With the same compactness tools as in the above lines, we conclude from estimates (5.25gi, (15.25hb . 
and (EBUl that there exists G L r (0, T; W rl - r (fi\r c )) n BV([0, T]; W iy (fl\T c )*) s.t. 

A d and A -± d in i r (0, T; W hr (n\T c )) for 1 < r < -, (6.7a) 
tf T , # T ->-tf ini r (0,T;W rl - £ ' r (fi\r c ))forallee (0, 1]. (6.7b) 

The latter convergence yields that 1 T , § T i? in L r (0, T; L 15 / 7 - e (ft)) for all e G (0, 8/7]. Taking into 
account estimate (5.25g i in L°°(0, T; L 1 (f2)) and arguing by interpolation, we conclude 

T , -dr-^d in L 5/3 - £ (Q) for all e G (0, §] . (6.7c) 

Notice that, under condition ( 14.7b on 9q, ( 15.7b and d6.7cb imply the strict positivity of *&. Moreover, Helly's 
selection principle and the a priori bound for (d T ) T in L°°(0, T; L 1 (J7)) yield 

■& T (t) A fi(t) in M(fi) for all i G [0, T]. (6.7d) 

Combining ( I6.7cb with ( 14.2b . it is immediate to deduce 

9(tf T ) 6(tf) inL 2 (Q). (6.8) 

Furthermore, it follows from d6\7bl the trace theorem d ^ [tf] : W 1 '^ (Ci\T c ) L 10 / 7 - £ (r c ) for all 
e G (0, 3/7], and gill, that 

[6(? T )J [0(0)] ini™(0,T;i"( 10 / 7 - e )(r c )) VeG (0, f ] . (6.9) 

Exploiting (16. lib . ( 16.2b . and (I6.1kb . which in particular yields 

lim/ z T (t)a ([u r (t)])dS= [ z(t)a ([u{t)})dS for all t G [0, T], 

we conclude by lower semicontinuity arguments that 

<f>(u(t),z(t)) < lirninf $ T (ur(t),Zr(t)) for all i G [0,T1. (6.10) 

T— >0 

Step 1: passage to the limit in the momentum equation. At first, we take the limit as r — » of the 

discrete momentum equation ( 15.12b with smooth test functions v G C°°((3; R 3 ), fulfilling [u] ^ on 
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S c . We approximate them with discrete approximations {v^} such that Ji^'J >: on S c and the related 
piecewise constant and linear interpolants fulfill, as r — s- 0, 



v T ^v mW 1 ' 1 {Q,T;L 2 {n; 
v T ^v mL 2 (0,T;W 2 - 2 (n\T c ;R 3 )), } 

ll Ve («r)|| L7( Q ;R 3x3x3 ) < C. 

In order to pass to the limit in the first integral term on the left-hand side of (15.121 . we use ( 16. lab . (16. lbb . 
and (16.8b . combined with ( 16.1 lb . The regularizing 7-terms in ( 15.121 ) vanish in the limit due to ( 16.111 ). Notice 
that d6.1kb and the continuity of a' (cf. (14.1 H i) imply a^Jur]]) «o([[u]) in L°°(E C ). Therefore, taking 
(16. 2b into account, we deduce that z^agdlvj]) — x zq;q([uJ) in L°°(S C ). This convergence and d6.1 lb 
allow us to pass to the limit in the second integral term on the left-hand side of d5.12b . To take the limit 
of the third and fourth terms (in the case g > 0), we use d6.1cb , d6.1db , and (16. let , as well as the first of 
(16. 1 lb . Combining the latter with d5.2b we also take the limit of the first term on the right-hand side of 
(15. 12b . The convergence of the other two integrals ensues from d5.9b , d6.1ab . and (16. lib . Thus, we have 
proved that the triple (u, z, 1?) fulfills equation (13.12b with smooth test functions. With a density argument, 
we conclude (EHJ) with test functions v G i 2 (0, T; ^(fi\r c ; R 3 )) n If 1,1 (0, T; L 2 (0; I 3 )). For later 
convenience, let us observe that, in the case K(x) is a linear subspace for almost all x € T c (cf. (14.5b ). 
taking test functions v = u + Xv with v any admissible test function satisfying [u]] >z on S c and A an 
arbitrary real number, we obtain ( 13. 12b in the form 

gu{T)-v(T)dx + / (De(it)+Ce(it)-Be(??)):e(w) + (HVe(w)+GVe(u)):Ve(i>) - gu-vdxdt 



+ I a' ([u}y[v}dSdt= I guQ-v{Q)dx + [ F-vdxdt - [ f-vdSdt (6.12) 
for any v G L 2 (0, T; W^ 2 (ft\r c ; K 3 )) n W 1 ' 1 ^, T; L 2 (ft; M 3 )). 

Step 2: passage to the limit in the semistability condition. We consider a subset N c (0, T) of full 
measure such that for all t € IN" the approximate stability condition (15 . 1 4b holds for the (countably many) 
considered r's. Then we fix £ e N and z G L°°(r c ). We may suppose without loss of generality that 

%(u(t),z—z(t)) < +00, hence 

z(x)<z(t,x) for a. a. x £ T c ■ (6.13) 
Then, we construct the following recovery sequence 

!z{x) 

z T (t,x)—, r where z(t, x) > 0, 
TK >z(t,x) V ' (6.14) 

where z(t, x) = 0. 

Now, using ( 16. 1 3b and ( 16.2b one immediately sees that 

z T {t, ■) <z T (t,-) a.e. inr c , z T {t) ^- z m L°° {T c ). (6.15) 

Plugging z T in d5.14b and using ( 16. lkb . (16.2b . and ( 16. 15b . we find 

< lim U T {v n .{t),Z T {t))+^{v n .{t),Z T {t) -Z^{t))-^ T {Ur{t),Zr{t)) 

lim / (a ([u r (i)])-a 1 ([u r (i)]))(S T (i)-z r (t))d5 



C 



(a ([u(t)j)- ai (lu(t)]))(z(t)-z(t))dS 



c 



= <Z>(u{t),z(t))+3l(u(t), z(t)-z{t))-$(u(t),z(t)) . (6.16) 

Step 3: passage to the limit in the mechanical and total energy inequalities. Using ( 16. lab . d6.1eb . ( 16.3b . 
and ( 16. 10b . we pass to the limit on the left-hand side of the discrete mechanical energy inequality d5.16b by 
weak lower semicontinuity. To take the limit of the right-hand side, we employ d5.2b (which in particular 
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yields 3> r (u 0jT , z ) —> 3>(u , Zo)), the weak convergence (16.1 al l and the strong convergence (16.8b , which 
give 

9(^ r )B:e(?V) 6(tf)B:e(u) weakly in i^Q). (6.17) 

We pass to the limit in the two remaining terms by (15.9ab - (15.9bb . Hence, the triple (u, z, «?) complies for 
all t G [0, T] with 

T kin (w(f)) +$(u(T),z(f)) + / / ie(u(s)):e(i(s))+GVe(«(s)):Ve(u(s)) dzds 

Jo Jn 

+ Var3j(u,z; [0,i]) < T kin (u ) + $ (u 0) z ) 

+ [ ( f Q(#(s))M:e(u(s))dx+ [ F(s)-u(s)dx + [ f(s)-u(s) ds) ds . (6.18) 
Jo \Jn Jn Jr N J 

By the very same lower semicontinuity arguments (also using ( 16.7db and ( 15.1 Oi l), we also pass to the limit 
in the discrete total energy inequality ( 15.131 l. 

Step 4: mechanical energy equality. First of all, we observe that the following chain rule-type inequality 
holds for all t G [0, T] 

<t>(u(t),z(t)) -$(u ,z )+Vaix(u,z;[0,t}) > f (X,u) ds 

Jo 

for any A G L 2 (0, T; W 2,2 (fi\r c ; K 3 )*) with X(t) G d u $(u(t), z(t)) for a.a. t G (0,T), (6.19) 

where d u $ : W 2 ' 2 (ft\T c ; R 3 ) =1 W 2 - 2 (il\r c ; K 3 )* denotes the subdifferential w.r.t. u of the functional 
$ : W^ 2 ^ 2 (r2\r c ;R 3 ) x L°°(r c ) ~> R defined in < f3~8T >. Easy calculations show that the operator d u $ is 
given by 

Ae<9„$(u,z) if and only if 3£e d"3 K (u); Vu G W 2,2 (fi\r c ;R 3 ) : 

(X,v) = [ Ce(u):e{v)+U\7e(uy.\7e(v)dx+ [ za' ([uj )• [vj dS + (£,v), (6.20) 
Jn Jr c 

where, for notational convenience, we have introduced the functional 3k ■ W 2 ' 2 (£l\T c ; R 3 ) — > [0, +oo] 
defined by 3 K (u) = I K ( \u\ ) (cf. (gj)), and its subdifferential d3 K : W 2 ' 2 (fl\T c ; R 3 ) =4 W 2 ' 2 (0\r c ;R 3 )* 
In order to prove ( 16. 1 9b for a fixed selection X(t) G d u <i>(ii(i), z(i)), we exploit a technique, combining 
Riemann sums and the already proved semistability condition (13.14b . which is well-known in the analy- 
sis of rate-independent systems and dates back to 1PJ. The main difficulty here is to adapt such a trick 
to the case of a Stieltjes integral (cf. (16.25b ). and to do so we will mimick the argument in the proof 
of (37] Prop. 3]. For any n > 0, we take a suitable partition = t% < V£ < . . . < t^ n = T 
with maxj=i jv„(£™— < 1/". in such a way that the functions A n : [0,T] — > L°°(r c ) given 
by A n (t) := oi([w(*?_i)]) for i G (*?_!,*?] fulfill 

yU->Oi(H) inL°°(E c ) asn->oo. (6.21) 

The existence of such partitions follows from the fact that u : [0, T] — > W 2 ' 2 (tt\T c ; R 3 ) is continuous, 
since u G W 1,2 (0, T; W 2 ' 2 (il\T c ; R 3 )). Thus it is also uniformly continuous, and so is the mapping 
JwJ : [0, T] — > i°°(r c ). Then, we use that uniformly continuous mappings admit uniform approximation 
by piecewise constant interpolants. In fact, (16.21| i holds for all partitions of [0, T] whose fineness tends 
to 0, therefore we can choose our partition in such a way that the semistability (13.14b holds at all points 
{*" : i = 0, . . . , N n —1, n G N}. Hence, we write (13.14b at f£_ x tested by z = z(tf), thus obtaining 

sKtSU), *(*?_!)) < + a(«(t?_ 1 ),*(t?)-*(*SLi)) 

= *(«(*?),*(<?)) + / axCHiti)])^^)-^*"-!)^^- f' (X n (s),u(s))d s 

for any selection A„(t) G d u G>(u(t), z(tf)) for a.a. i G tp], « = 1, N n , where we have also used 

the chain rule for the convex functional u M- z(tf)), cf. 11481 Prop. XI.4. 11]. In particular, taking into 



24 



RICCARDA ROSSI AND TOMAS ROUBICEK 



account formula d6.20t for we choose 

X n (t) = A(t) - p n (t) withp„(t) G W^ 2 (n\T c ;R 3 )* given by 

(Pn(t),v) := / (z(t)-z(tf))a' ([u(t)])- {vj dS. (622) 
Summing for i = 1, . . . , N n , we obtain 

<S>(u(T),z(T))-<S>(u ,zo)+J2 _ A n \z\(dSdt)>J2 (X n (s),u(s)) ds. (6.23) 
Now, reproducing the calculations throughout ||35l Formulae (4.70)-(4.74)], it can be shown that 

liminfV / (X n (s),u(s))ds> (X(s),u(s)) ds. (6.24) 
On the other hand, it follows from d6.21b that 

lim V f' / .A n |z|(dSdi) = / ai([u] )z(dS,di) = Var K (u,z; [0,T]). (6.25) 

Indeed, i G C(E C )* can be extended to i°°(S c )* by the Hahn-Banach principle, and then tested by 
A n - ai({u\) G L°°(£ c ) which converges to zero by d6.21b . Combining ( l6.23M6.25b , we obtain (16.19b . 
In order to make d6.19b more explicit, we may observe that 

(£,u) ds = J K (u(t)) - J K (u(0)) = I K ([u(t)})- Ik{{u(0)]) = 

(6.26) 

foralH G L 2 (0,T;I^ 2 ' 2 (O\r c ;M 3 )*) such that i(s) G 9J K (u(s)) fora.a.s G (0,T), 

by the chain rule for the convex functional 3k (cf. again [48, Prop. XI.4.11]), and by d3.7al i and (13.11b . 
Therefore, in view of ( I6.19b - d6.26b , we conclude the following inequality for all t G [0, T] 

$(u(t), z{t)) - $(w , zq) + Vara(u, z; [0, t]) 



> 



o \ Jn 



Ce(w):e(w)+MVe(u):Ve(w)da;+ / zao([u])-[uJ dSJds. (6.27) 



In order to develop the test of (13.12b by ii, we need to distinguish the quasistatic case g = and the 
dynamical case g > 0. 

Case g > 0. First of all, let us observe that, under (14.5b . the qualification v G VK^ 2 (0, T; L 2 (fi; R d )) for 
the test functions in d3.12b might be relaxed to 

v G L 2 (0,T;^ 2 (r2\r c ;M 3 ))nW^ 1 ' 2 (0,T;^ 2 (fi\r c ;R d )*), (6.28) 

cf. notation (13.5b . Indeed, thanks to (13.11b and to the linearity of K(x) for almost all x G T c , the function u 
fulfilling ( |3~T2l is such that ii G L 2 (0, T; W^ 2 (^\r c ; R d ))- Note that d6T28T > is sufficient to give meaning 
to the term J Q ii-v dxdt, because the spaces L 2 (0, T; W^ 2 (fl\r c ; R d )*) and L 2 (0,T;W K ' 2 (n\T c ;R d )) 
are in duality. 

Now, a comparison in (I6TT21 yields that u G L 2 (0, T; W / ^' 2 (^\r c ; R 3 )*). Therefore, (@3]i ensues, and 
ii is an admissible test function for the momentum balance inclusion (13.12b , since it fulfills (16.28b . Then, 
upon proceeding with such a test we conclude for all t G [0, T] that 



, u(t)\ 2 dx+ [ [ De(u):e(u)+GVe(M):Ve(w)dxds 
n Jo J n 



2 



+ Ce(u):e(u)+EX7e(u)'Ve(u)dxds+ / / za' ( [u] )• [it] dSds 
Jn Jn Jo Jr c 

= | / |M | 2 da;+ f ( [ 0(D)W.e(u) dx + f F-udx+ f f-uds) ds. (6.29) 
2 Jn Jq V Jn Jn Jr N J 



ADHESIVE CONTACT DELAMINATING AT MIXED MODE, ITS THERMODYNAMICS AND ANALYSIS 



25 



Combining (16.29b with ( 16.27b , we get the converse of inequality d6.18b , hence the desired mechanical 
energy equality ( I3.18l l ensues. 

Case g = 0: A comparison in (13.12b with g = shows that the functional 



e-.v^ (Be(u)+Ce(u)-B0(t?)):e(t;) + (GVe(w)+HVe(V)):Ve(u) dccdi 
Jq 

+ I za' (lu}y[vj dSdt - [ F-vdxdt- [ f-vdSdt (6.30) 
Js c JQ Je n 

is in L 2 (0, T; W 2 ' 2 (Cl\T c ; R 3 )*), and fulfills 

I K ([v])dt> f I K ([u})dt+ [ (£,v-u)dt. (6.31) 
Jo Jo 

Hence, £(t) € dJ K (u( t)) for almost all t £ (0, T). Thus, (ET26t yields /„' (£, u) ds = for all i G [0, T], 
which is just relation d6.29b with g = 0. Again, we combine the latter with (16.27b . and conclude the 
mechanical energy equality (13.181 ). 

Step 5: passage to the limit in the enthalpy equation. First of all, we observe the following chain of 
inequalities for all t € [0, T]: 

Vax%(u,z;[0,t]) + [ [ Be(u):e(u)+GWe(u)':We(u) dxds 
Jo Jn 

<liminf/ / d ({vJ , z T ) dSds + f f Be(^):e(^)+GVe(^):Ve(iV) dxds 
Jo Jr c Jo Jn 

< hmsupT kin (u 0iT ) + $ T (w ,r, Zq) ~ Tkin(wrW) - $ T (tv(t), Z T (t)) 

+ f ( f e(5 T )B:e(ir) + / F r -?Vd.T+ / 7r'«r ds) ds 
Jo \is! Jn Jr N / 

< T kin (u ) + $(uo,2 ) - T kin (w(i)) - #(u(t),z(t)) 

6(t?)B:e(ti) +F-udx+ f f-udS)ds 
Jr N / 

= Diss s (u,z; [0,*]) + / / De(u):e(u)+GVe(w):Ve(u) dzds. (6.32) 
Jo Jn 

Indeed, the first inequality ensues from ( 16. 1 ab and d6.3l ). the second one from the discrete mechanical energy 
inequality ( 15.16b . the third one from (EU, ( 16. lei ). ( 16.101 ). ( 16.17b . and from d5.9a| )— ( f579bb , cf. also Step 3. 
Finally, the last equality ensues from (13.18b proved in Step 4. Thus, all of the above inequalities turn out to 
hold as equalities. By a standard liminf/limsup argument, we find in particular 

ne(ur)-e{ur) -> De(u):e(u) and GVe(u T ):Ve(u r ) -> GVe(u):Ve(u) strongly in L^Q). 

Combining these convergences with ( 16.17b we pass to the limit in the first term on the right-hand side of 
( 15.15b . To take the limit of the second right-hand-side term, we observe that 

lim J ai([u r ] )z T vdSdt = lim J (ai([tEr] ) — ai([ur] ty&rVdSdt 

+ lim [ a 1 ([v T ])z T vdSdt = + j oi([u])u«(d9dt) (6.33) 

for any i> e C(E C ), and in particular for v = — ^ ; here we used respectively ( 16. 5b , ( 16.6b , and ( |5.25fb . 

Then, we pass to the limit in the left-hand side of ( 15.15b by exploiting ( I6.1bl ). ( I6.7al ). d6.7bl ). d6.7db . d6.9b , 
as well as properties d4.1cb for % and d4.1el ) for r\, and by arguing in the very same way as in the proof 
of BTl Thm. 5.1], to which we refer for all details. 

At the end, employing (15. 10b . we take the limit of the last term on the right-hand side of d5.15b . thus 
finding that the triple (u, z, 1?) fulfils the weak formulation d3.15b of the enthalpy equation. 
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Step 6: total energy identity. We test the weak formulation ( 13.151 ) of the enthalpy equation by 1 and add 
it to the mechanical energy equality. This gives the total energy balance (13.13b . □ 
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